On the vertex corrections in antiferromagnetic spin fluctuation theories 
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We argue that recent calculations by Amin and Stamp 
(PRL 77, 301 (1996)) overestimate the strength of the vertex 
corrections in the spin-fermion model for cuprates. We clar- 
ify the physical origin of the apparent discrepancy between 
their results and earlier calculations. We also comment on 
the relative sign of the vertex correction. 



In a recent publication Q Amin and Stamp computed 
the leading order correction to the spin-fermion vertex in 
the nearly antiferromagnetic Fermi-liquid model. They 
considered a region near optimal doping in which the 
Fermi surface is large, and precursors of the spin-density- 
wave state have not been formed yet. They argued that 
the leading vertex correction is rather large (Ag/g ~ 1) 
for realistic values of the parameters of the electron dis- 
persion and the spin fluctuation spectrum. This result 
is in clear disagreement with previous calculations 0-@] 
which reported a much smaller amplitude of the cor- 
rection. Whether vertex corrections are large or not is 
relevant to the validity of the Eliashberg-type calcula- 
tions of the superconducting transition temperature in 
cuprates ||. 

In this letter, we clarify the origin of the discrepancy 
between Ref. [1] and earlier results and show that Ref. 
[1] strongly overestimates the amplitude of the vertex 
correction near optimal doping. 

The spin-fermion model describes fermions coupled to 
spin fluctuations by 
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Here g is the coupling constant, and <7j are the Pauli 
matrices. It is assumed that the Fermi-liquid descrip- 
tion is valid, i.e., near the Fermi surface, the electronic 
Green's function behaves as G(k,uJ m ) = Z/(iuj m — e^), 
where tk = £k — jU, and Z < 1 is a positive con- 
stant. The dispersion near the Fermi surface is given 
by £fc = — 2i(cos k x + cos k y ) — At' cos k x cos k y . Spin fluc- 
tuations are described by a dynamical spin susceptibility 
which is assumed to be strongly peaked near the anti- 
ferromagnetic momentum Q — (tt, tt), and to behave at 
low energies as x{l->^m) = Xq/0- + £ 2 <7 2 + ^mlAv)- 
Here q 



and oj 



of 



= ? — Qi £ is the magnetic correlation length, 
oc £~ 2 is a typical spin fluctuation frequency 



the problem due to the proximity to antiferromagnetism. 
Amin and Stamp computed the vertex corrections for 
the points on the Fermi surface which are connected by 
Q (the so-called "hot spots"), and also for k along the 
Brillouin zone diagonal. They argued that the relative 
correction is larger for the hot spots - a result we do 
not dispute. For the vertex correction at the hot spot 
they obtained assuming that fermions are free particles 
Ag/g = -{g 2 XQ^ s f/4:ir 3 n 2 )I(k h ) where I(k h ) is a di- 
mensionless quantity. It turns out that I{kh) is rather 
large which implies a large amplitude of the correction. 
Amin and Stamp claim that the origin of the discrep- 
ancy between their and earlier results lies in their 
more accurate numerical evaluation of I(kh)- We dis- 
agree with their explanation and will argue that the rea- 
son for their discrepancy with earlier results in fact has 
a physical rather than a numerical origin. 

We independently computed the vertex corrections an- 
alytically and numerically along the same lines as de- 
scribed above. The analytical expression for Ag/g is 
given by 
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J{2-f v £) = £uj s f/v <C 1, and v and 



which is much smaller than any other energy scales in 



where 5 = 

define the fermionic dispersion around a hot spot: efe = 
v(k — kh) cos(0) and eu+Q — v(k — kh) cos(</> + </>q). The 
angle </>o is the angle between the normals to the Fermi 
surface at the hot spots (see the Fig.p. 

We emphasize that Eq. (||) contains a term which 
logarithmically depends on 6 and a term independent 
of 6. Previous analytical calculations of the vertex cor- 
rection restriced with the | log (5 1 term only. We will 
see that for parameters relevant to cuprates both terms 
nearly equally contribute to the vertex renormalization. 

We first compare our analytical (Eq.(|^)) and numeri- 
cal ||] results with the ones obtained by Amin and Stamp 
for the same set of parameters, namely t = Q.25eV, t' — 
-0.45i, \fx\ = 1.46i, uj sf = UmeV, £ = 2.5a, g = OMeV 
and Z = I. They found Ag/g w -0.7 while Eq.(§) 
yields Ag/g w 0.55, and numerical computations || yield 
Ag/g ~ 0.4 (for this set of parameters <&o = 1-78 and 
5 = 0.27). We see that apart from the sign difference 
which we discuss later, our results and the ones by Amin 
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FIG. 1. The graphical representation of the angle 0o be- 
tween the normals to the Fermi surface at hot spots (dashed 
lines). For clarification we omitted the parts of the Fermi 
surface in the second and fourth quadrant. 



and Stamp are in close agreement which implies that a 
more accurate evaluation of I[kh) cannot be the main 
reason for the discrepancy mentioned above. The actual 
reason for the discrepancy with earlier work lies in the 
fact that Amin and Stamp considered u> s f as an input 
parameter, whose value can be inferred from NMR ex- 
periments, and simultaneously treated fermions as free 
particles with Z = 1 (!). We will show that these two 
assumptions are likely to be incompatible. 

The point is that near optimal doping, the key source 
for the spin damping is the interaction with fermions (the 
damping due to this interaction is much stronger than 
the one due to a direct spin-spin exchange). Since the 
damping term in the spin susceptibility is related to the 
imaginary part of the particle-hole bubble at momen- 
tum transfer Q and low frequencies, the fermions which 
contribute to uj s f are located in the vicinity of the hot 
spots. One can then compute the damping within the 
spin-fermion model of Eq.(|l|) which, we remind, is valid 
only near the Fermi surface where fermionic excitations 
are coherent. The result is 
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If one substitutes the same values for the parameters as 
above, one obtains ui s f ~ l.OGmeV which is more than 
ten times less than ui s f = lAmeV used in the above cal- 
culations of Ag/g. A way to restore the experimentally 
inferred value of u> s f is to assume that the fermions are 
not free particles, i.e., Z < 1. Specifically, one needs 
Z = 0.28 to obtain the correct value of ut s f. Substituting 
Z = 0.28 into the formula for Ag/g (Eq.(||)), we imme- 
diately obtain Ag/g « +0.04. Note in passing that a 
small value for Z is consistent with photoemission data 
Q which show that even at optimal doping, the quasi- 
particle peak is rather broad, and there is a substantial 
spectral weight in the incoherent part of the spectrum. 

Alternatively, one can consider fermions with arbitrary 
Z, compute uj s f using Eq.(0), and substitute the result 



into Eq.(g). Doing this, we obtain neglecting terms of 
0(6") 

Ag | sin 0o | f" log[sin(0/2)] sin(0/2) 
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This is the way the vertex correction was obtained in [Q . 
For the set of parameters chosen, we indeed recover the 
same result. Similar calculations for 0o = T have been 
performed in Ref. Q which also yielded a very small value 
of Ag/g. 

We now address two technical points. The first one is 
the sign of the vertex correction. We have already men- 
tioned that the sign obtained in earlier calculations J§-|J] 
is opposite to the one obtained by Amin and Stamp. A 
simple way to check the sign of the correction is to con- 
sider the limit where both the coupling constant and the 
chemical potential are much larger than the fermionic 
bandwidth, and w g/2. In this limit the electronic 
structure develops precursors of the spin-density-wave 
state, as two of us have recently demonstrated explic- 
itly Accordingly, the renormalized spin-fermion ver- 
tex should be much smaller than the bare one due to a 
"near" Ward identity (our argument here parallels the 
one recently displayed by Schrieffer &). Meanwhile, a 
simple examination of Eq.(6) in Ref. jl| shows that at 
large /x, the vertex correction is positive, i.e., Ag/g — > +1 
rather than tending to —1 which is necessary to obtain 
the physically motivated strong reduction of the vertex. 

The second point concerns the computations near the 
edge of the antiferromagnetic instability, when the corre- 
lation length is very large. It follows from Eq. (^) that as 
6 oc ujgfl; oc £ x tends to zero, Ag diverges logarithmically. 
In this situation, higher-order corrections are indeed rel- 
evant. It has been shown in |^,||] that the logarithms sum 
up to a power law, and the full vertex takes the form 
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As we already discussed, however, (3 is negative and nu- 
merically quite small. We thank D. Pines and P.C.E. 
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